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Results on the Gluon Propagator in Lattice Covariant Gauges 

L. Giusti a , M. L. Paciello b , S. Petrarca* b c , C. Rebbi a , B. Taglienti b 

a Boston University - Department of Physics, 590 Commonwealth Avenue, Boston MA 02215 USA 
b INFN, Sezione di Roma 1, P.le A. Moro 2, 1-00185 Roma, Italy 
c Dipartimento di Fisica, Universita di Roma "La Sapienza" 

We illustrate preliminary results on the gluon propagator computed in generic covariant lattice gauges in the 
quenched approximation with the Wilson action. We have applied a recently proposed procedure to fix a generic 
covariant gauge on the lattice. We make a comparison among results obtained at two different values of the gauge 
parameter and in the standard Landau's gauge. 



We present preliminary data on the gluon prop- 
agator in generalized covariant gauges. The mo- 
tivation for considering a generalized gauge-fixing 
procedure stems, just as in the continuum, from 
the necessity of studying the gauge dependence of 
the Green's functions. In the covariant gauges it 
is possible to change the gauge fixing varying the 
value of the gauge parameter. In the continuum, 
the generic covariant gauges satisfy the following 
condition: 



(1) 



where A\(x) are matrices belonging to the Lie al- 
gebra of the SU(3) group and A is the gauge pa- 
rameter which controls the width of the gaussian 
distribution of the A's. 

The corresponding gauge fixing procedure on 
the lattice is implicitly defined by the following 
expression for the expectation value of a gauge 
dependent operator O: 

(O) = ' 
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(2) 



where S(U) is the Wilson lattice gauge invariant 
action and G\ is the gauge transformation that 
enforces the gauge condition. 

After having generated a set of A according to 
a Gaussian distribution, the gauge fixing proce- 
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dure follows very closely the widely used Lan- 
dau's gauge fixing: the gauge transformation 
G(x) which rotate the links U^x) into the gauge- 
fixed configuration is obtained by the numerical 
minimization of a functional -F[G], chosen in such 
a way that its minimum corresponds to a gauge 
configuration satisfying the gauge condition. The 
functional form we adopt in the continuum reads: 



F[G] = J d i xTr [(d^A 



A) 2 



(3) 



The discretization of this expression requires a 
particular care (driven discretization) in order to 
obtain a functional form suitable to be minimized 
by the standard gauge- fixing algorithms. More 
details on this covariant gauge-fixing method can 
be found in 0, @. 

The stationary points of the above functional 
correspond to the following gauge condition in the 
continuum: 



A,&,(0 M A?f*-A) = O 



(4) 



It should be noted that, while for A = cq. (jj) re- 
duces to the Landau gauge condition (apart from 
the presence of possible spurious solutions, i.e. so- 
lutions to d^A Gx — A = p with D v d v p — which, 
however, we did not encounter in our study), the 
same is not true of the lattice construction, where 
the equations that follow from the minimization 
of the functional with A = differ from the lattice 
Landau gauge condition by terms 0(a). Thus the 
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comparison of the results we will obtain for A = 
with those obtained in the Landau gauge will be 
a meaningful check on our procedure. 

We gauge-fixed 221 SU(3) configurations with 
V = 16 3 x 32 and (3 = 6 taking the gauge pa- 
rameter A = 0, 8 where A = has been con- 
sidered in order to make a comparison with the 
Landau gauge. The gauge configurations were re- 
trieved from the repository at the "Gauge Con- 
nection" (http://qcd.nersc.gov/). The calcula- 
tion has been performed on Boston University's 
Origin2000 and took approximately 19000 CPU 
hours. We imposed a gauge fixing quality factor 
9 < 10~ 6 . With these gauge-fixed configurations 
we computed the correlator of gluon propagator 

D^ix-y) = <A^{x)A u {y) > (5) 

where we adopted for A the standard expression 

~U,(x)~U^x) 



2iag 



(G) 



- Traceless 

Discussions on the effect of different lattice defini- 
tions of the gluon field in the context of numerical 
calculations can be found in refs. ||] and The 
general properties of lattice operators were stud- 
ied in ref. ||. 

As usual in gluon propagator calculations j(| 
we study the quantity 

1 



(AoAo)(t) 



V 2 ^ 



Tr(4)(x, iL4 (y,0)> . (7) 



which is displayed in Fig. 111. This correlator 
should be a constant in the case of Landau's 
gauge-fixing and its behaviour is indicative of the 
quality of the gauge fixing procedure. Figure |l| 
shows that even in the case A = 8 the above cor- 
relator is almost a costant. The case A = is 
equivalent to the Landau's one with a different 
value of the constant. 

In Fig. | we exhibit the behaviour of the prop- 
agator 

(AA>(t) = ^YlJ2 Tr (M*,t)My,o)) (8) 

i x.y 

as a function of t. 

The results of this simulation show clearly that 
the gauge dependence of the gluon propagator 
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Figure 1. Behaviour of {AoAo)(t) as a function of 
time for a set of 221 thermalized SU (3) configurations 
at /3 = 6.0 with a volume V x T = 16 3 x 32. The data 
have been slightly displaced in t to help the eye, the 
errors have been obtained with the jacknife method. 



(eq. (H)) on A cannot be simply understood as a 
rescaling of the data, as one might have argued of 
the data on a smaller lattice, reported in the very 
preliminar analysis of ref. 0]. We also observe 
that the data for the Landau gauge fixing and for 
A = are compatible within errors. We proceed 
by showing the gluon propagator in Fourier space. 
The transformed field A^ (k) is given by 



i^^e-^-^A^x) 



(9) 



where = j^-n^ with L M = 16 or 32 for spatial 
and temporal indices respectively, are integers. 
Following suggestions from lattice perturbation 
theory Q, in order to get rid of a well-known 
lattice artifact we express all the quantities as a 
function of the variable q 



?At = 2sin(-^) 
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Figure 2. Gluon propagator eq. (JSJ) as a function of t 
for a simulation with 221 thermalized configurations 
at f3 = 6 and V — 16 3 x 32. The apparent asymmetry 
with respect to the midpoint of the lattice is only due 
to the slight shift of two sets of points along the time 
axis. 



Figure 3. Transverse part of the gluon propaga- 
tor q 2 DT(q) in covariant gauges as a function of q. 
The two sets of data refer to A = and A = 8, 221 
thermalized SU(3) configurations at f3 — 6.0 with a 
volume V x T = 16 3 x 32. 



rather than k. This kinematic correction is 
needed to reduce the anisotropy of the data when 
displayed as function of k. The Fourier transform 
of the gluon propagator in eq. ^ is: 
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D„u(q) =y< A^AS (q) > 



(10) 



We analyzed the data in terms of the conventional 
distinction of transverse and longitudinal parts 
Dt and D^: 



)D T {q) 



q^qv Dbiq) 



(11) 



q , r r 

In Fig. [| we present the behaviour of the trans- 
verse part (times q 2 ) of the gluon propagator in 
covariant gauges, for the two values of the gauge 
parameter A — 0, equivalent to the Landau's 
gauge, and A = 8, as a function of q = y/q 2 . The 
data have been averaged over the Z3 group and 



only values corresponding to ko = ||- (0 -r- 12) and 
ki = 7j|(0-j-6) are reported. We see that at large 
q 2 the data at the two values of the gauge param- 
eter tend to coincide, while the gluon propagator 
at low q 2 exhibits a clear gauge dependence. The 
data with A = 8 show a lower peak than those 
with A = 0, but the peak position stays fixed. In 
Fig. U we illustrate the behaviour of the longitu- 
dinal part as a function of q. 

The data at A = vanish as one expects for 
Landau's gauge. At A = 8 the situation is differ- 
ent and a signal different to zero is observed with 
a clear increase for increasing q. This differs from 
the predictions of pertubation theory, where the 
longitudinal part is a constant proportional to the 
value of the gauge parameter. 

In summary, our gauge fixing procedure can be 
implemented for medium to large volume simula- 
tion; the gauge fixing obtained with our algorithm 
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Figure 4. Longitudinal part of the gluon propagator 
as a function of q at A = and A = 8. 



behaves as one expects in the case of A = 0, re- 
producing the features of Landau's gauge. The 
data of the gluon propagator are new and show 
a clear gauge dependendence; the transverse and 
longitudinal part will be analyzed in more details 
in a future publication. 
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